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Abstract A brief review of relativistic effects in few-body systems, of theoretical ap- 
proaches, recent developments and applications is given. Manifestations of relativistic 
effects in the binding energies, in the electromagnetic form factors and in three-body 
observables are demonstrated. The three-body forces of relativistic origin are also dis- 
cussed. 
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1 Introduction 

Light nuclei give a typical example of a few-body weakly bound system. Their binding 
energies B are of the order of 0.1% from their masses M. This however does not mean 
that the relativistic effects in light nuclei are also so small: they are much larger than 
the ratio B/M. The reason is that, in contrast to the hydrogen atom (for instance), 
the small nuclear binding energy is a results of cancellation of much more significant 
kinetic and potential energies. This can be illustrated as follows. 

Let us consider a system of two non-relativistic particles interacting by the Yukawa 
potential: V non — re i, (r) = — ^- exp(— fir) that in the momentum space corresponds to 
the kernel: 

Vnon — rel.iO) 9 '. o- (-0 

This non-relativistic kernel is a limiting case of the following relativistic one-boson 
exchange kernel: 

, . — 47TQ — 4:Tva . . 

V re i.{q) = — 2 ~ = — 2 2 ( 2 ) 
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Fig. 1 Left: Binding energy B vs. coupling constant a, found via Shrodinger equation (NR), 
Bethe-Salpeter (BS) and the light-front (LFD) one, for exchange mass fi = 1 in the kernels JTJ 
and J2]l. Right: the same as at left, but for different exchange masses fi (here the BS and LFD 
results are indistinguishable from each other). The results are from [4]. 



V re i turns into V non _ re ; when go — ^ 0. The kernel V non _ re i enters the Shrodinger 
equation. The kernel V re i_ enters the relativistic Bethe-Salpeter (BS) equation |1 . An- 
other popular relativistic approach is light-front dynamics (LFD). The corresponding 
equation - light front (LF) equation, see for review 2,3 , - contains the kernel V re i , 
which analytical form differs from but its non-relativistic limit also coincides with 
Q. All these three equations - Shrodinger equation, BS and LF ones - were solved, 
for spinless particles in the J — state, in Ref. [4]. The results in the limit of ex- 
tremely small binding energy B — > are shown in fig. [T] Left panel corresponds to 
heavy exchange mass fj, = 1 (in the units of the constituent mass m = 1). We see 
that the relativistic calculations BS and LF are very close to each other and, at the 
same time, they strongly differ from the non-relativistic result NR. This shows that 
the relativistic effects can be important even at small binding energy. The curves at 
the right panel demonstrate that when the exchange mass decreases, the difference be- 
tween relativistic and non-relativistic results decreases too. That is, the true relativistic 
system is a system not only with very small binding energy but also with interaction 
resulted from exchange by zero mass. When the exchange mass is not small (of the 
order of the constituent mass m), then the particles are in a very narrow potential 
well of the radius r ~ 1/m. Then their momenta are comparable with their masses 
k ~ 1/r ~ m, the kinetic energy (and the system at all) is relativistic and, hence, 
the small binding energy is a result of cancellation of large (positive) kinetic and large 
(negative) potential energies. That's why a system with small binding energy may be 
still relativistic. Similar situation is realized in nuclei, since the exchange mesons like 
U) and p are enough heavy ~ /i p ~ 0.8 m). 

In the opposite case of strongly bound systems, there is an example [5] that even 
a system with extremely large binding energy, such that its total mass tends to zero, 
may be dominated (by 90%) by a few lowest Fock sectors, containing two, three and 
four particles. It still remains to be a few-body system (though, highly relativistic). 

Even in the case of mainly non-relativistic system (average momentum is very 
small), its impulse distribution contains a relativistic tail. This tail may be very small, 
but it completely determines the e.m. form factor at large momentum transfer. The 
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form factor of this system is also very small but it should be calculated in a relativistic 
approach. 

All that requires development of appropriate approaches to relativistic description 
of few-body systems. Brief review of these approaches is given in the next section. 

2 Relativistic descriptions 

In non-relativistic quantum mechanics the wave function is an eigenvector of Hamilto- 
nian: Hifj = Etjj. Dynamics is introduced by the adding to free Hamiltonian an 
interaction term H mt = V: 



In relativistic case, the relativistic covariance is guaranteed if the wave function 
(or the state vector |p)) is forming a representation of the Poincare group. The lat- 
ter is determined by ten generators P^ , J M „ which satisfy the following commutation 
relations: 



The Hamiltonian H = Pq is now only one of generators. Similarly to non-relativistic 
case, dynamics is introduced by the adding to free Poincare generators an interaction 
terms in a way which keeps the commutation relations @ unchanged. This is not simple 
but a solvable task. It can be realized in the framework of two different approaches: (i) 
relativistic quantum mechanics with fixed number of particles and (ii) field theory. 

In relativistic quantum mechanics the Poincare generators are the functions of fixed 
number (say, two or three) of the particles momenta. The interaction is a phenomeno- 
logical one, it is fitted to describe e.g. the two-body phase shifts. Then one can make 
predictions: to calculate for instance the e.m. form factors or three-body observables. 
For good reviews of this approach see J6][T] . 

In field theory the Poincare generators are derived from Lagrangian by a well- 
known formulas, given almost in any textbook. If Lagrangian is not free (contains 
interaction), then the interaction appears also in the generators. The state vector \p), 
on which the generators act, can be decomposed in the basis of free fields (similarly to 
the Fourier decomposition in plane waves of the non-relativistic wave function). This 
basis is represented as an (infinite) set of Fock components with increasing numbers of 
particles. In practice, this decomposition is truncated (the desired number of particles 
is fixed by hand). After that the approach becomes approximate. 

So, in practice, two approaches - (i) and (ii) - differ from each other by the point 
where one makes this truncation: (i) either from the very beginning, with further 
phenomenological construction the generators; (ii) or after finding the generators by 
the field theory recepees. In the latter case, the kernel is motivated by a field-theoretical 
Lagrangian. In its turn, this field-theoretical interaction is mainly restricted by the one- 
boson-exchange. 

One should also distinguish three forms of relativistic dynamics, proposed by Dirac 
[B], which exist in both approaches. Namely: (a) instant form; (b) front form; (c) 
point form. They differ from each other by the ways of introducing the interaction in 
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generators. In the instant form the time component Pq of the four-momentum operator 
Pn contains interaction, whereas the spatial components Pj (j = x,y,z) are free. The 
interaction enters also in the components Jqj of the operators J^„. The components 
Jy are free. In the front form (LFD) the interaction enters the component P— = 
Po — P z , whereas the components P+ — Pq + P z , P x , P y are free. The operator J^j, 
is constructed correspondingly. The components of this operator which transform the 
LF plane t + z — const into itself are free. In the point form all the components of P^ 
contain interaction, whereas the operator J^u is free. 

In its turn, LFD is developed in two forms: ordinary LFD with the LF plane 
t + z — const (see for review [3]) and explicitly covariant LFD 12 with the LF plane 
given by the invariant equation u) ■ x — uj^t — oj • x = const, where ui = (uq,u:) is a 
four-vector with cj 2 = 0. The main advantage of this latter formulation is in the fact 
that the dependence of the state vector on the LF orientation is given explicitly, in 
terms of the four- vector uj (see e.g. the LF deuteron wave function Q below). In the 
particular case u = (1, 0, 0, —1) we come back to the ordinary formulation of LFD. 

The operator of e.m. current used to calculate e.m. form factors, like any four- 
vector operator, has the following commutation relation with J Kp : 



If J Kp contains interaction (in l.h.-side of ©), then r.h.-side of i.e. j^, also must 
contain interaction. This means that in interacting system the exact e.m. current cannot 
be free (except for the point form of dynamics). 

Another series of (field-theoretical) relativistic approached deals not with the state 
vector \p) itself, but it is based on the BS amplitude 1. defined as: 



tp(x) is the Heisenberg field operator and (0| is the vacuum state. In the momentum 
space: 

$ = $(ki,k 2 ,p) = $(k,p), p 2 = A/ 2 , fc 2 /m 2 , fcl/m 2 and fc = (fe 1 -fc 2 )/2. 

The BS equation for <& is singular, that complicates its numerical solution. To avoid 
singularities, one can transform this equation in the Euclidean space. Corresponding 
solution provides the binding energies. However, to calculate e.m. form factors, we 
should know the BS amplitude in Minkowski space. The Wick rotation in terms of the 
relative momentum k - the argument of the BS amplitude <P(k,p) - is not valid in the 
form factor integral over k. 

The methods to solve BS equation in Minkowski space were recently developed 
first for the spinless particles 9 and then for two fermions [10] , They are based on the 
Nakanishi integral representation [llj : 



This integral determines a singular BS amplitude <&(k,p). However, the Nakanishi 
weight function 9(7,2) is not singular. Substituting BS amplitude in the form (J6]l in 
the BS equation, one can derive and solve numerically equation for 9(7, z). Then again 
using © with known 51(7,2), one can express the observables, like e.m. form factors, 
through g("/,z) analytically and then compute them numerically. Another method to 
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solve the BS in Minkowski space is based on the separable approximation of the kernel 
(see [12] and references therein). 

There are also a few reductions of the four-dimensional BS amplitude to a three- 
dimensional form (still in the Minkowski space). In this direction, the approach pro- 
posed in [13] (covariant spectator theory) is most advanced and well developed. In the 
covariant spectator theory, the NN potential was fitted and applied to the deuteron 
and three-body problems as well as to the e.m. form factors [14, 15 . 

The theoretical activity in studying the relativistic few-body systems flourishes in 
all the forms of dynamics and in all the approaches listed above. 

3 Some applications 

It is clearly demonstrated (see e.g. |14j ) that the non-relativistic calculations of the ed 
elastic cross section do not describe the data at Q 2 > 1 GeV 2 /c 2 . One needs to perform 
the calculations with true relativistic deuteron wave function. In the non-relativistic 
case the latter is determined by two spin components: S- and D-waves. In relativistic 
approaches the number of components depends on approach. In the spectator theory 
[13j . there are four components. There are six components in covariant LFD [TS]. The 
deuteron BS amplitude is determined by eight components (see e.g. [2]). 

As an example, we mention the calculation carried out in the framework of explicitly 
covariant LFD. In this approach, the relativistic deuteron wave function has the form 

EM- 

V(k,n) = / X -L<7 + Ja|( 3k fr; g) - a) + / 3 |(3n(n. a) - a) (7) 

+/ 4 ^(3k(n • a) + 3n(k ■ a) - 2(k ■ n)a) + h\J\\[* x n] + /e|f [[k x n] x a], 

where n = and er are the Pauli matrices. The vector n just provides the explicit 

dependence of this wave function on the LF orientation. The six components /i_6 
were calculated in ]T7]. Only three of them dominate: /i,/2 (which turn into the S- 
and D-waves in non-relativistic limit) and f§, whereas /3,/4,/6 are negligible. The 
corresponding deuteron e.m. form factors were calculated in [18] . The results of this 
calculation are in good coincidence with the appeared later experimental data [19] . 
We do not give this comparison here. A detailed review can be found in IT. One 
can conclude that the relativistic effects in a two-body system are taken into account 
satisfactory. 

On the contrary, there are still the problems in the theoretical descriptions of the 
three-body systems. Though the problems with binding energy of tritium (underbind- 
ing) can be removed by incorporating the three-body forces, there are some deviations 
in description of the elastic pd scattering. They are seen in fig. [2] taken from [201121] . 

There is also a discrepancy in the analyzing power A v in the pd elastic scattering 
(see left panel in fig. [3] taken from [25]). Right panel [23 shows the nd breakup cross 
section: nd — > (nn)p in a particular kinematics corresponding to so-called symmet- 
ric space-star configuration. In both cases, there are considerable deviations between 
different versions of the theoretical calculations and experimental data. One hopes to 
resolve these contradictions, properly taking into account relativistic effects as well 
as three-body forces. We will see below that the three-body forces can be partially 
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Fig. 2 Left: Cross section of elastic nd scattering vs. cm. scattering angle. Short dashed line 
- non-relativistic calculation. Solid line - non-relativistic + 3-body forces. The figure is taken 
from 20 . Right: The same as at left for other energies. Solid line - non-relativistic calculation. 
Dashed line - relativistic one. The figure is taken from 1211 . 




Fig. 3 Left: Analyzing power in elastic pd scattering vs. cm. scattering angle. Dotted line — 
non-relativistic calculation. Solid line — relativistic one. The figure is taken from 1221 . Right: 
The cross section of the nd breakup: nd — > (nn)p. The figure is taken from 1231 . 

induced by relativity. The importance of relativistic effects in exclusive pd breakup 
scattering at intermediate energies was demonstrated in 24, 25], where the relativistic 
Faddeev equation was solved without employing a partial wave decomposition. The 
relativistic effects improve agreement with experimental data. The magnitude of these 
effects depends on configuration in the final state. Some success in describing A y was 
achieved in [25]. 



4 Relativity in three-body systems 

The binding energy of two-body system interacting by a potential described by the 
potential well Uq with radius ro tends to constant when Uo — > oo, ro — > but Uor^ = 
const. On the contrary, for this interaction, the binding energy of three-body system 











3 , 





















Fig. 4 Left: Three-body binding energy B vs. coupling constant a, found via Shrodinger 
equation (long-dashed), Bethe-Salpeter (solid) and the light-front (short-deshed), for exchange 
masses /i = 0.01 and fi = 0.5. The results are from 1301 . Right: The graphs contributing in the 
three-body forces of relativistic origin. 



tends to — oo. This is a well-known property of non-relativistic three-body system which 
is called the Thomas collapse |27| . 

It turned out that the relativity results in an effective repulsion: for given two-body 
mass M2, the three-body mass M3 is finite 28,29 . This drastic change of three-body 
binding energy shows that the influence of relativity on three-body system may be 
stronger than on the two-body one. 

However, for enough strong interaction, corresponding to large two-body binding 
energy, such that M2 < M c ~ 1.43 m, the three-body mass becomes negative. In this 
domain of M2, a physical solution for the three-body system disappears. In the non- 
relativistic scale, the binding energy equal to the total mass of constituents, is almost 
infinity. Therefore the case, when M3, though being finite, approaches to zero, is a 
relativistic counterpart of Thomas collapse. 

The relativistic three-body equations - BS and LF ones - have been also solved, for 
spinless particles, not only for zero range interaction, but also in more realistic case of 
one-boson exchange [50]. The corresponding two-body solution for binding energy [3] 
was discussed above and is shown in fig. [1] The three-body binding energy vs. coupling 
constant a is shown in fig. [4] (left panel). In contrast to the two-body results (see 
fig. [T]), the BS and LF calculations do not coincide, but considerably differ from each 
other. They both also differ from the non-relativistic result (like in the two-body case). 
However, for the two-body system, when the exchange mass /j, tends to zero, the BS 
and LF calculations (which are very close to each other) tend to the non-relativistic 
result. In three-body system this is not the case. The reason of these deviations is 
the three-body forces generated by relativity. Corresponding graphs, containing two 
mesons in flight (first considered in [31]), are shown in the right panel of fig. [4] They 
are automatically included in the three-body BS equation. However, they should be 
added explicitly in the kernel of the LF equation. After taking them into account 30 , 
we find good coincidence between the BS and LF results (see fig. [5J. This explicitly 
demonstrates that in a three-body system (a) relativistic effects and three-body forces 
appear together; (b) both may be important. Notice that the role of three-body forces 
may be different in different relativistic approaches. Thus, the relativistic three-body 
forces are not generated as a correction to the three-body spectator equation [131115] . 
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Fig. 5 Three-body bound state mass squared JV/| vs. coupling constant a for exchange masses 
fi = 0.01 (left), fj, = 0.5 and fi = 1.5 (right). The units are set by the constituent mass: m = 1. 
The figures are taken from 1301 . 



However, they should be still incorporated as a relativistic correction to the Schrodinger 
equation. 

The relativity is not the only source of the three-body forces. There exist other 
sources (e.g., the intermediate isobar creation) which may generate "intrinsic" three- 
body forces. One should include all that in the analysis of the discrepancies in three- 
body reactions discussed above. 



5 Conclusions 

We conclude that relativistic effects in nuclei can be important in spite of small bind- 
ing energy. At high momenta they clearly manifest themselves and are necessary to 
describe the deuteron e.m. form factors. At the same time, there is still a discrepancy 
in three-body observables which might be a result of less clarity in understanding the 
corresponding relativistic effects, the relativistic NN kernel and the three-body forces. 

Relativistic few-body physics remains to be a field of very intensive and fruitful 
researches. 



References 

1. E.E. Salpeter, H.A. Bethe, Phys. Rev. 84, 1232 (1951). 

2. J. Carbonell, B. Desplanques, V.A. Karmanov, J.-F. Mathiot, Phys. Rep. 300, 215 (1998). 

3. S.J. Brodsky, H.-C. Pauli, S.S. Pinsky, Phys. Rep. 301, 299 (1998). 

4. M. Mangin-Brinet, J. Carbonell, Phys. Lett. B 474 237 (2000). 

5. Dae Sung Hwang, V.A. Karmanov, Nucl.Phys. B 696) 413 (2004). 

6. F. Coester, Prog, in Part, and Nucl. Phy s. 29, 1 (1992). 

7. W.N. Polyzou, J. Lab lectures. larXiv:0908.14"4Tl (2009). 

8. P.A.M. Dirac, Rev. Mod. Phys. 21, 392 (1949). 

9. V.A. Karmanov, J. Carbonell, Eur. Phys. J. A 27, 1 (2006). 

10. J. Carbonell, V.A. Karmanov, Eur. Phys. J. A 46, 387 (2010). 

11. N. Nakanishi, Phys. Rev. 130, 1230 (1963); Graph Theory and Feynman Integrals (Gordon 
and Breach, New York, 1971). 

12. S.G. Bonda renko, V.V. Burov , W.-Y. Pauchy Hwang, E. P. Rogochaya, Nucl. Phys. A 832, 
233 (2010); arXiv:0810.4470 (nucl-th'): larXlvT l002.0487 (nucl-th). 

13. F. Gross, Phys. Rev. 186, 1448 (1969); Phys. Rev. D 10, 223 (1974); C 26, 2203 (1982). 

14. R. Gilman, F. Gross, J. Phys. G 28 R37 (2002). 

15. S.A. Pinto, A. Stadler, F. Gross, Phys. Rev. C 79, 054006 (2009); 81, 014007 (2010); 



9 



16. V.A. Karmanov, Nucl. Phys. A 362, 331 (1980). 

17. J. Carbonell, V.A. Karmanov, Nucl. Phys. A 581 625 (1995). 

18. J. Carbonell, V.A. Karmanov, Eur. Phys. J. A 6, 9 (1999). 

19. D. Abbott et al., Phys. Rev. Lett. 84, 5053 (2000). 

20. H. Witala, W. Glockle, D. Ruber, J. Golak, H. Kamada, Phys. Rev. Lett. 81, 1183 (1998). 

21. H. Witala, J. Golak, W. Glockle, H. Kamada, Phys. Rev. C 71, 054001 (2005). 

22. H. Witala et al., Phys. Rev. C 77, 03400 4 (2008). 

23. H. Witala, W. Glockle, arXiv:0909.0495 (nucl-th). 

24. T. Lin, Ch. Elster, W.N. Polyzou, W. Glockle, Phys. Lett. B 660, 345 (2008). 

25. T. Lin, Ch. Elster, W.N. Polyzou, H. Witala, W. Glockle, Phys. Rev. C 78, 024002 (2008). 

26. F. Gross, A. Stadler, Phys. Lett. B 668, 163 (2008). 

27. L.H. Thomas, Phys. Rev. 47, 903 (1935). 

28. T. Frederico, Phys. Lett. B 282, 409 (1992). 

29. J. Carbonell, V.A. Karmanov, Phys. Rev. C 67, 037001 (2003). 

30. V.A. Karmanov, P. Maris, Few-Body Syst. 46, 95 (2009). 

31. Shin-Nan Yang, Phys. Rev. C 10, 2067 (1974). 



